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ABSTRACT 

Quantifying  air-sea  exchanges  of  enthalpy  and  momentum  is  important  for  understanding  and  skillfully 
predicting  tropical  cyclone  intensity,  but  the  magnitude  of  the  corresponding  wind  speed-dependent  bulk 
exchange  coefficients  is  largely  unknown  at  major  hurricane  wind  speeds  greater  than  50  m  s-1.  Since 
direct  turbulent  flux  measurements  in  these  conditions  are  extremely  difficult,  the  momentum  and  en¬ 
thalpy  fluxes  were  deduced  via  absolute  angular  momentum  and  total  energy  budgets.  An  error  analysis  of 
the  methodology  was  performed  to  quantify  and  mitigate  potentially  significant  uncertainties  resulting 
from  unresolved  budget  terms  and  observational  errors.  An  analysis  of  six  missions  from  the  2003  Coupled 
Boundary  Layers  Air-Sea  Transfer  (CBLAST)  field  program  in  major  hurricanes  Fabian  and  Isabel  was 
conducted  using  a  new  variational  technique.  The  analysis  indicates  a  near-surface  mean  drag  coefficient 
CD  of  2.4  X  10  3  with  a  46%  standard  deviation  and  a  mean  enthalpy  coefficient  CE  of  1.0  X  10  3  with 
a  40%  standard  deviation  for  wind  speeds  between  52  and  72  m  s-1.  These  are  the  first  known  estimates  of 
CK  and  the  ratio  of  enthalpy  to  drag  coefficient  CK/CD  in  major  hurricanes.  The  results  suggest  that  there  is 
no  significant  change  in  the  magnitude  of  the  bulk  exchange  coefficients  estimated  at  minimal  hurricane 
wind  speeds,  and  that  the  ratio  CKICD  does  not  significantly  increase  for  wind  speeds  greater  than 
50  m  s_1. 


1.  Introduction 

Heat,  moisture,  and  momentum  exchanges  at  the  air- 
sea  interface  are  primary  processes  in  tropical  cyclone 
(TC)  intensification  and  maintenance  (Malkus  and  Riehl 
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I960;  Emanuel  1986, 1995).  Wind  speed-dependent  bulk 
aerodynamic  formulas  often  are  used  to  represent  the 
turbulent  fluxes  associated  with  air-sea  interactions, 
relying  on  drag  CD ,  moisture  CE ,  and  enthalpy  CK  ex¬ 
change  coefficients.  However,  as  Ooyama  (1969)  stated, 
“Unfortunately,  there  is  little  information  onC£  under 
hurricane  conditions,  other  than  the  semispeculative 
guess  that  the  exchange  coefficients  of  latent  heat,  sen¬ 
sible  heat,  and  momentum  are  probably  of  the  same 
magnitude.”  Nearly  40  yr  after  that  pioneering  study, 
enthalpy  and  momentum  exchange  coefficients  are 
still  largely  unknown  at  major  hurricane  wind  speeds 
(>50  m  s~\  equivalent  to  category  3  and  higher  on  the 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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Saffir-Simpson  scale).  Since  the  ocean  surface  charac¬ 
teristics  change  significantly  with  increasing  wind  speed 
(Black  et  al.  1986),  a  simple  extrapolation  of  the  co¬ 
efficient  magnitudes  derived  at  low  wind  speeds  to 
50  m  s-1  and  beyond  is  not  necessarily  justified.  The 
current  research  aims  to  improve  our  understanding  of 
air-sea  interaction  in  major  hurricanes,  given  the  high 
impact  that  their  damaging  winds  and  storm  surge  can 
have  on  coastal  populations,  global  economics,  maritime 
activities,  and  naval  operations. 

Emanuel  (1986)  derived  an  expression  for  the  po¬ 
tential  intensity  (hereafter  E-PI)  of  a  TC  that  depends 
on  local  environmental  parameters  and  is  proportional 
to  the  ratio  of  the  bulk  enthalpy  and  momentum  ex¬ 
change  coefficients  CK/CD.  Montgomery  et  al.  (2006) 
and  Bell  and  Montgomery  (2008,  hereafter  BM08) 
showed  that  a  30  m  s”1  range  of  E-PI  estimates  for 
category  5  Hurricane  Isabel  (2003)  could  be  obtained  by 
varying  the  magnitude  of  this  ratio  from  0.5  to  1.5.  More 
refined  estimates  of  the  bulk  exchange  coefficients 
would  clearly  help  narrow  the  uncertainty  in  E-PI  pre¬ 
dictions.  The  importance  of  understanding  of  air-sea 
exchange  at  high  wind  speeds  is  not  limited  to  the  the¬ 
oretical  domain,  however.  E-PI  estimates  are  used  in 
statistical  TC  intensity  forecasts  (Emanuel  et  al.  2004), 
and  numerically  simulated  intensity  has  shown  sensi¬ 
tivity  to  prescribed  surface  exchange  coefficients  in 
a  variety  of  simplified  model  frameworks  (Ooyama 
1969;  Rosenthal  1971;  Rotunno  and  Emanuel  1987; 
Emanuel  1995).  Three-dimensional,  full  physics  models 
have  also  shown  sensitivity  to  the  parameterization  of 
surface  fluxes  in  hurricane  intensity  forecasts  (Braun 
and  Tao  2000;  Bao  et  al.  2002;  Davis  et  al.  2008).  It  must 
be  noted  that  parameterizing  the  complexity  of  the  air- 
sea  interaction  at  these  wind  speeds  with  10-m  bulk 
exchange  coefficients  is  likely  an  oversimplification,  but 
accurate  estimates  are  still  important  given  the  estab¬ 
lished  practical  reliance  on  these  coefficients.  Improve¬ 
ments  in  numerical  weather  models’  representation  of 
air-sea  interaction  in  major  tropical  cyclones  would 
likely  contribute  to  improved  intensity  forecasts  issued 
by  the  Joint  Typhoon  Warning  Center  and  the  National 
Hurricane  Center  (NHC)  (Rappaport  et  al.  2009). 

Surface  fluxes  at  major  hurricane  wind  speeds  are  very 
difficult  to  measure  in  the  tropical  cyclone  boundary 
layer  (TCBL).  Given  that  the  spatial  and  temporal  oc¬ 
currence  of  these  wind  speeds  is  very  small  compared  to 
the  ocean  basin  size  and  seasonal  time  scales  in  which 
hurricanes  occur,  the  probability  of  a  fixed  sensor  en¬ 
countering  them  is  very  low.  Instrumented  ships  pru¬ 
dently  avoid  these  conditions  in  general,  and  operational 
ocean  sensors  (such  as  buoys)  are  not  designed  to 
withstand  the  extreme  conditions  found  in  the  open 


ocean  at  the  base  of  a  major  hurricane  eyewall.  There¬ 
fore,  laboratory  experiments  and  observations  from 
research  aircraft  deployed  in  tropical  cyclones  are 
currently  the  most  viable  choices  for  obtaining  infor¬ 
mation  about  surface  fluxes  at  high  wind  speeds.  How¬ 
ever,  low-altitude  flying  in  the  turbulent  boundary  layer 
at  the  base  of  an  eyewall  is  also  very  hazardous,  which 
requires  that  remote  sensing  or  expendables  (i.e., 
dropwindsondes)  be  used  to  obtain  measurements  in 
this  region. 

Global  positioning  system  dropwindsonde  (hereafter 
simply  dropsonde)  profiles  have  been  used  to  estimate 
flux  profiles  and  CD  at  the  highest  wind  speeds  yet 
(Powell  et  al.  2003;  Vickery  et  al.  2009).  Powell  et  al. 
(2003)  provided  the  first  indications  of  a  slight  decrease 
in  CD  from  maximum  values  around  2.5  X  10  3  at  30-40 
m  s"1.  Dropsonde  profiles  analyzed  by  Vickery  et  al. 
(2009)  indicated  a  continued  slight  decrease  in  CD  at 
wind  speeds  up  to  60  m  s_1.  Additional  evidence  that 
CD  does  not  increase  above  35  m  s_1  was  reported  by 
Donelan  et  al.  (2004)  using  laboratory  tank  measure¬ 
ments.  The  eddy  correlation  method  was  used  for 
wind  speeds  up  to  26  m  s-1,  and  a  momentum  budget 
retrieval  was  used  for  wind  speeds  from  20-50  m  s_1. 
These  tank  measurements  showed  a  “saturation”  of  CD 
around  hurricane-force  wind  speed  (33  m  s_1)  and 
suggested  a  limiting  aerodynamic  roughness  of  the  sur¬ 
face  waves  above  these  speeds.  Indirect  retrievals  of  CD 
in  hurricanes  that  are  consistent  with  the  above  studies 
have  also  been  conducted  using  ocean  measurements 
(Shay  and  Jacob  2006;  Jarosz  et  al.  2007). 

The  uncertainty  in  the  magnitude  of  CD  and  CK  in 
the  TCBL  was  one  of  the  issues  investigated  as  part 
of  the  Coupled  Boundary  Layers  Air-Sea  Transfer 
(CBLAST)  field  campaign  (Black  et  al.  2007).  Because 
of  the  safety  hazards  mentioned  previously,  instru¬ 
mented  P-3  aircraft  were  only  flown  in  the  clear  air 
boundary  layer  between  rainbands,  but  fortunately  fa¬ 
vorable  conditions  for  turbulence  observations  were 
achieved  in  two  major  hurricanes,  Fabian  and  Isabel 
(2003).  French  et  al.  (2007)  and  Drennan  et  al.  (2007) 
reported  the  first  open-ocean  eddy  correlation  mea¬ 
surements  in  the  TCBL  at  strong  tropical  storm  force 
wind  speeds  (20-30  m  s_1).  Estimates  of  CD  and  CE 
from  these  studies  indicated  no  discernible  dependence 
on  wind  speed  in  the  range  measured.  Although  the 
French  et  al.  (2007)  results  slightly  differ  from  those 
reported  by  Powell  et  al.  (2003)  and  Donelan  et  al. 
(2004),  general  agreement  exists  that  an  extrapolation  of 
an  increasing  CD  from  low  wind  speed  formulas  (e.g., 
Large  and  Pond  1981)  is  not  justified.  Recent  research 
suggest  that  CK  is  nearly  independent  of  wind  speed 
from  20  to  30  m  s_1  using  the  CBLAST  measurements 
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(Zhang  et  al.  2008)  and  from  13  to  40  m  s”1  using  lab¬ 
oratory  measurements  (Haus  et  al.  2010;  Jeong  et  al. 
2012).  There  have  been  no  previous  estimates  of  the 
magnitude  of  CK  at  major  hurricane  wind  speeds  to  the 
authors’  knowledge. 

Emanuel  (2003)  developed  a  similarity  theory  for 
gradient  wind-based  exchange  coefficients  at  very  high 
wind  speeds  and  suggested,  based  on  an  asymptotic 
analysis,  that  both  exchange  coefficients  should  approach 
constant  values  in  the  limit  of  high  wind  speed.  He  also 
presented  a  mechanistic  argument  that  suggests  that  the 
ratio  of  enthalpy  to  momentum  coefficients  should  de¬ 
crease  with  increasing  ocean  temperature.  But  neither 
of  these  approaches  predicts  specific  values  of  the  ex¬ 
change  coefficients.  The  scaling  theory  of  Emanuel  (2003) 
is  based  on  the  gradient  wind  and  neglects  the  explicit 
role  of  the  (unbalanced)  boundary  layer  dynamics  and 
their  contribution  to  the  total  wind  speed.  The  total  wind 
speed  is  the  intensity  measure  of  most  interest  to  hur¬ 
ricane  forecasters. 

An  indirect  approach  to  retrieving  momentum  ex¬ 
change  using  an  absolute  angular  momentum  budget 
was  originally  proposed  over  50  yr  ago  (Palmen  and 
Riehl  1957).  In  their  formulation,  the  surface  stress  is 
obtained  by  measured  transports  of  absolute  angular 
momentum  in  an  axisymmetric  cylindrical  coordinate 
system.  If  the  flux  of  angular  momentum  is  known  at  the 
top  and  sides  of  a  prescribed  control  volume,  the  flux  at 
the  air-sea  interface  can  be  obtained  via  residual,  as¬ 
suming  that  the  system  is  in  a  steady  condition.  Addi¬ 
tional  studies  used  angular  momentum  budgets  to 
deduce  CD  with  improved  datasets  (Miller  1962,  1964; 
Hawkins  and  Rubsam  1968;  Hawkins  and  Imbembo 
1976).  The  quality  and  density  of  hurricane  observations 
have  improved  significantly  since  that  time,  as  well  as 
the  analysis  techniques  used  to  composite  the  data.  The 
central  focus  of  our  research  was  the  application  of 
a  similar  methodology  to  retrieve  surface  fluxes  in  major 
hurricanes  using  recent  observations.  In  this  study,  we 
utilize  budgets  of  absolute  angular  momentum  and 
total  energy  to  retrieve  CD  and  CK  in  major  hurricanes 
Fabian  and  Isabel. 

In  practical  application,  some  of  the  momentum  and 
energy  budget  terms  are  very  difficult  to  calculate  with 
observational  data,  and  uncertainties  must  be  acknowl¬ 
edged.  To  gain  insight  into  the  relative  importance  of  the 
various  sources  of  error  in  the  budgets,  the  sensitivity  to 
errors  in  unresolved  budget  terms,  sea  surface  tem¬ 
perature  (SST),  center  location,  gridding  method,  and 
the  size  of  budget  control  volumes  were  examined 
systematically  using  numerical  model  data  where  the 
surface  fluxes  were  known.  The  results  of  the  sensitivity 
tests  were  then  utilized  to  determine  quantitative 


estimates  of  the  uncertainty  of  the  magnitude  of  the 
retrieved  bulk  exchange  coefficients. 

A  high-resolution  dataset  collected  in  Hurricanes 
Fabian  and  Isabel  in  2003  as  part  of  the  CBLAST 
experiment  (Black  et  al.  2007)  was  used  to  apply  an 
energy  and  momentum  budget  method  to  real  tropical 
cyclones.  Both  TCs  achieved  category  4  intensity,  and 
the  data  collected  represent  a  significant  advance  in 
the  quantity  and  quality  of  observations  in  major 
hurricanes.  To  deduce  momentum  and  enthalpy  fluxes 
from  these  observations  using  the  current  methodology, 
the  data  are  gridded  in  an  axisymmetric  coordinate 
system  moving  with  the  tropical  cyclone.  Montgomery 
et  al.  (2006)  and  BM08  used  a  simple  objective  analysis 
scheme  (Barnes  1973)  to  derive  a  gridded  axisymmetric 
structure  for  Hurricane  Isabel  for  use  with  the  E-PI 
calculations.  A  new  variational  analysis  procedure  called 
Spline  Analysis  at  Mesoscale  Utilizing  Radar  and  Air¬ 
craft  Instrumentation  (SAMURAI)  was  used  here  to 
improve  the  derived  axisymmetric  structure  used  in 
the  budget  calculations. 

A  review  of  the  equations  used  in  the  budget  calcu¬ 
lations  and  summary  of  the  uncertainties  are  presented 
in  section  2,  with  more  detailed  descriptions  of  the  error 
analysis  and  SAMURAI  technique  provided  in  the  ap¬ 
pendices  for  interested  readers.  The  results  of  the 
CBLAST  data  analysis  are  presented  in  section  3,  fol¬ 
lowed  by  the  derived  air-sea  fluxes  and  bulk  exchange 
coefficients  in  section  4.  The  concluding  section  sum¬ 
marizes  the  key  findings  of  this  study. 

2.  Methodology 

a.  Conservation  of  absolute  angular  momentum 

To  determine  the  bulk  momentum  exchange  coeffi¬ 
cient,  we  utilize  the  conservation  of  azimuthally  averaged 
absolute  angular  momentum  in  cylindrical  coordinates, 
defined  as  M(r,z)  =  rv  +  (l/2)/r2,  where  r  denotes  the 
radius  from  the  storm  center,  v  denotes  the  tangential 
(swirling)  wind  velocity, /denotes  the  Coriolis  parameter, 
and  the  overbar  represents  both  an  azimuthal  and  tem¬ 
poral  average.  For  the  current  analysis,  we  hold  /to  be  a 
constant  in  the  azimuthally  averaged  coordinate  sys¬ 
tem.  The  equation  for  conservation  of  M  in  flux  form  is 
given  by 

a(pM)  |  djpruM)  |  d(pwM)  d(j\e)  r^[ze  _  Q 
dt  rdr  dz  rdr  dz 

(i) 

where  the  overbars  have  been  dropped  for  clarity.  The 
Reynolds  stresses  are  defined  as  rre  =  -pu'v'  and 
Tze  =  -pw'v\  where  the  subscripts  represent  the  wind 
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directions  of  the  covariances  ( r  is  the  radial  direction,  6 
is  the  tangential  direction,  and  z  is  the  vertical  direction), 
primes  represent  the  departures  from  the  azimuthal  and 
time  averages,  p  denotes  the  density,  u  denotes  the 


radial  wind  velocity,  and  w  denotes  the  vertical  velocity. 
Integrating  over  a  control  volume  from  z i  to  z2  and  rx  to 
r2  yields  an  integrated  flux  form  over  a  control  volume: 


From  Eq.  (2),  the  radially  integrated  surface  stress 
may  be  calculated  from  measurements  on  the  sides  and 
top  of  an  arbitrary  control  volume.  A  schematic  of  a 
hypothetical  control  volume  is  shown  in  Fig.  1,  which 
illustrates  the  approximate  locations  of  zi,  z 2,  C,  and  r2 
for  an  idealized  hurricane  flow.  We  assume  that  the 


vertical  stress  at  z  1  can  be  represented  by  a  bulk  formula 
rze=  —  pw'v’  =  pCc>|u/J  u,  where  |u^|  denotes  the  hori¬ 
zontal  wind  speed,  and  neglect  the  radial  variation  of  CD 
over  the  control  volume  to  yield  the  drag  coefficient  in 
terms  of  mean  quantities,  with  unresolved  terms  grou¬ 
ped  into  a  residual  labeled  R: 


CD  [  r 1 


CZ2  fZ2  (r2  (r2  \  /  fr2  . 

ij  [puM]  |  ^  dz  ~  r2  [puM]  |  ^  dz  +  J  [pwM]\ zr  dr  -  J  [pwM]\ zr  dr  +  R  j  /  [p\uh\v]\ zr  dr , 


(3) 


where  the  unresolved  residual  is 


Alternatively,  one  can  obtain  the  average  surface 
stress  in  the  control  volume  by  dividing  Eq.  (2)  by 
JJ2  r2  dr ,  or  obtain  the  average  friction  velocity  u*  by 
dividing  Eq.  (2)  by  -  fr*pr2  dr  and  taking  the  square 
root.  Note  that  z  1  does  not  necessarily  have  to  be  at 
10-m  altitude  but  could  be  anywhere  in  the  surface 
layer  assuming  the  fluxes  are  nearly  constant  in 
that  layer.  If  z  1  is  higher  than  10  m,  then  the  de¬ 
nominator  in  Eq.  (3)  must  be  calculated  using  the 
10-m  density  and  winds  to  obtain  the  surface  drag  co¬ 
efficient. 

b.  Conservation  of  total  energy 

To  determine  the  bulk  moist  enthalpy  exchange  co¬ 
efficient,  we  utilize  the  conservation  of  total  energy, 
given  in  flux  form  in  cylindrical  coordinates  by 


d(p£)  ,  djpmE)  ,  d{pwE)  JU 

dt  rdr  dz  dt  ^  \2  J 

+  kV2T  +  QR,  (5) 

where  the  total  energy  E(r,0,z,t)  =  cpT  +  Lq  + 
(1/2) (u2  +  v2  +  w2)  +  gz  is  the  sum  of  the  internal 
(sensible  and  latent  heat)  energy  and  mechanical  (ki¬ 
netic  and  potential)  energy  of  the  mean  flow.  Here  cp 
denotes  the  specific  heat  at  constant  pressure,  T  de¬ 
notes  the  temperature,  L  denotes  the  latent  heat  of 
vaporization,  q  denotes  the  water  vapor  mixing  ratio,  g 
denotes  gravity,  p,  denotes  the  molecular  viscosity,  k 
denotes  the  thermal  conductivity,  and  QR  denotes  dia- 
batic  heat  exchange  due  to  radiative  transfer.  Equation 
(5)  is  a  form  of  Bernoulli’s  equation  (Gill  1982).  Because 
the  thermodynamic  equation  was  expressed  in  terms  of 
the  enthalpy  instead  of  the  internal  energy  (e.g.,  cpT  = 
cvT  +  p/p),  a  local  time  derivative  of  p  appears  on  the 
right-hand  side  that  is  associated  with  acoustic  waves  in 
the  atmosphere.  Here,  we  assume  that  the  local  time 
derivative  of  p,  diffusive  component  of  friction,  thermal 
conductivity,  and  radiation  are  small  and  can  be  ne¬ 
glected  for  this  application.  Reynolds  averaging  Eq.  (5) 
over  azimuth  and  time  and  integrating  over  the  control 
volume  yields 
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\-Fzk\  \Zlrdr  = 


[Z2  - 

r2[puE  +  Frk  +  ue  +  u'e  -  wjrz  -  vTrg\\ r  dz 

L  2 

fZ2  - 

rx [puE  +  Frk  +  ue  +  u'e  -  wrrz  -  vt.J |  r  dz 

[pwE  +F,k  +  we  +  w'e  —  ut  —  vt  0] \  rdr 


[pwE  +  we  +  w'e  —  uttz  —  vrzd 


r  dr  dz, 


(6) 


where  all  variables  represent  azimuthal  and  temporal 
means,  with  the  bars  dropped  except  over  third-order 
covariance  terms  for  clarity.  The  turbulent  kinetic 
energy  is  given  by  e  =  (1/2 )p(u'2  +  v'2  +  w'2),  and  the 
eddy  enthalpy  fluxes  are  represented  by  the  F  terms,  where 
the  subscripts  indicate  the  respective  covariances  with  an 
oppositely  signed  convention  to  the  Reynolds  stresses. 
Note  that  dissipative  heating  terms  do  not  appear  in  the 
total  energy  balance  equation  (6)  since  they  represent 
a  conversion  between  mechanical  and  heat  energy,  nor  do 
terms  appear  involving  the  conversion  of  mean  kinetic 
energy  to  turbulent  kinetic  energy  (e.g.,  Lindzen  1990, 
90-92).  However,  flux  gradients  involving  the  in¬ 
teraction  of  turbulent  momentum  fluxes  and  shearing 
flow  remain.  For  lack  of  an  established  terminology, 
these  terms  are  referred  to  as  the  “shear  flux”  and 


represent  a  transport  of  energy  through  the  bound¬ 
aries  of  the  control  volume  by  the  interaction  of  turbu¬ 
lence  and  the  mean  flow.  For  the  lower  boundary, 
^shear  =  JJ2  [urrz  +  utz0]  represents  a  loss  of  kinetic 
energy  to  the  ocean  through  the  interaction  of  the  sur¬ 
face  momentum  fluxes  and  mean  flow.  The  shear  flux 
term  is  not  included  in  the  residual  since  it  may  be  resolved 
indirectly  via  the  derived  surface  momentum  fluxes. 

We  assume  that  the  vertical  moist  enthalpy  flux  at  z i 
can  be  represented  by  a  bulk  formula  Fzk  =  FzT  + 
Fzq  =  pCpW’T'  +  pLw'q'  =  Cxp\uh\  (k*  -  k),  where  k*  is 
the  saturation  moist  enthalpy  at  the  sea  surface  (k*  = 
cpTsst  +  Lq*)  and  q*  is  the  saturation  mixing  ratio  at 
the  surface.  Neglecting  the  radial  variation  of  CK  over 
the  control  volume  yields  the  bulk  enthalpy  exchange 
coefficient: 


C  K 


CZ2 

r2[puE]  | 

\\pwEMZx 


dz  ~ 


rdr  +  FSHE  AR  +  R 


\  [puE]  |  r  dz  +  J  [pwE]  \zr  dr 

J  \p\^h\(k*  -  k)]\ zr dr . 


(7) 


The  unresolved  residual  is 


[Z2 

R=  h[F, 


rk  +  ue  +  u'e-WTrz-vTr6 


zk  +  we  +  w'e  -  urn  -  vtz6 


fZ2  _ 

II  r2  llz~\_  rl  \-Frk  +  ue  +  u'e  ~  WTrz  ~  ”T  J I  rt  (lz 

'z2  p  \d(pE  +  e ) 


]|zrdr-J  [we  +  w'e]\  z  r  dr  +  J  J 


dt 


r  dr  dz.  (8) 


Alternatively,  one  can  obtain  the  average  enthalpy 
flux  by  dividing  Eq.  (6)  by  the  radially  integrated  ra¬ 
dius  similarly  to  the  average  momentum  flux.  Like¬ 
wise,  if  zi  is  higher  than  10  m,  then  the  denominator  in 
Eq.  (7)  must  be  calculated  using  the  10-m  density,  wind, 
and  enthalpy  to  obtain  the  surface  enthalpy  exchange 
coefficient. 


c.  Error  analysis 

Accuracy  in  the  retrieval  of  the  bulk  exchange  co¬ 
efficients  requires  (i)  minimal  errors  in  the  specification 
of  the  axisymmetric  kinematic  fields,  (ii)  strict  adher¬ 
ence  to  mass  continuity,  and  (iii)  a  small  magnitude  of 
the  unresolved  residual  terms.  Accurate  retrieval  of  CK 
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Fig.  1.  Schematic  illustrating  hypothetical  control  volume  (black 
dashed  line)  used  for  the  budget  methodology.  A  simplified  sec¬ 
ondary  circulation  (gray  streamlines)  and  a  contour  surrounding  the 
region  of  maximum  wind  (gray  oval  labeled  “Vmax”)  are  shown  to 
indicate  that  the  control  volume  encompasses  the  eyewall  region. 

also  requires  (iv)  minimal  errors  in  the  axisymmetric 
thermodynamic  fields  and  (v)  accurate  measurements 
of  the  SST.  To  test  the  accuracy  of  the  methodology, 
comprehensive  error  analyses  were  undertaken  using 
idealized  numerical  simulations.  A  summary  of  the  er¬ 
rors  identified  by  the  sensitivity  tests  is  presented  here, 
and  described  in  more  detail  for  interested  readers  in 
appendix  A.  Table  1  shows  the  magnitude  of  potential 
errors  that  could  affect  the  solution,  steps  that  were 
taken  to  mitigate  the  error,  and  the  estimated  uncer¬ 
tainties  remaining  after  mitigation.  If  left  unconstrained, 
the  unmitigated  errors  listed  can  easily  overwhelm  the 
magnitudes  of  the  retrieved  exchange  coefficients,  but 
they  are  small  enough  after  mitigation  to  allow  mean¬ 
ingful  interpretation  of  the  results.  The  largest  source  of 
potential  error  was  estimated  to  be  numerical  inaccu¬ 
racies  in  the  satisfaction  of  the  mass  continuity  equation. 
Strict  adherence  to  mass  continuity  was  enforced  by 
using  a  transverse  streamfunction  as  a  control  variable  in 
the  SAMURAI  technique,  thereby  eliminating  this 
source  of  error  from  the  solution.  The  second  largest 
source  of  error  results  from  the  neglect  of  the  unresolved 
residual  terms  in  the  budget.  The  volume-integrated 


tendency  term  and  vertical  eddy  fluxes  at  the  top  of  the 
control  volume  were  identified  as  the  largest  of  these 
unresolved  terms.  The  error  from  neglecting  these  terms 
manifests  primarily  as  sensitivity  to  the  size  and  shape  of 
the  control  volume.  An  optimal  set  of  control  volume 
sizes  was  identified  to  minimize  the  error,  and  the  re¬ 
maining  uncertainty  is  estimated  at  ±50%. 

The  shear  flux  term  is  about  40% -60%  of  the  mag¬ 
nitude  of  the  surface  fluxes,  which  results  in  a  significant 
low  bias  if  neglected.  However,  the  low  bias  can  be 
corrected  by  estimating  the  shear  flux  using  the  mean 
wind  multiplied  by  the  surface  stress  derived  from  the 
momentum  budget.  The  bias  correction  is  exact  where 
the  mean  winds  and  CD  are  known  exactly,  but  errors  in 
the  mean  winds  and  surface  stress  using  real  data  add 
some  uncertainty.  Sensitivity  tests  indicate  that  a  ±50% 
error  in  CD  translates  to  a  ±20%  error  in  the  magnitude 
of  CK  from  the  shear  flux  term.  To  avoid  adding  too 
much  noise  to  the  CK  estimates,  a  mean  CD  derived  from 
all  six  missions  was  used  in  this  study  to  estimate  the 
surface  stress  used  in  evaluating  the  shear  flux  for  CK. 
We  therefore  estimate  a  CK  uncertainty  of  ±20%  due  to 
the  bias  correction. 

The  sensitivity  tests  also  identified  analysis  errors, 
SST  errors,  and  circulation  center  errors  as  areas  where 
mitigation  was  required  to  ensure  the  best  possible  re¬ 
sults  from  the  observational  dataset.  Analysis  errors 
were  minimized  by  using  a  variational  approach  designed 
specifically  for  this  study  that  is  described  in  the  next 
section.  SST  errors  were  minimized  by  using  airborne 
expendable  bathythermograph  (AXBT)  data  in  combi¬ 
nation  with  satellite-derived  temperatures  as  described  in 
section  3b.  The  use  of  high-resolution  radar-derived  cir¬ 
culation  centers  minimizes  the  error  in  the  cylindrical 
coordinate  transform  as  described  in  section  3c.  The  re¬ 
maining  uncertainties  for  each  of  the  error  sources  after 
mitigation  are  estimated  to  be  on  the  order  of  ±20%. 
Fortunately,  these  errors  are  independent  and  random,  and 
there  is  no  reason  to  expect  that  they  would  be  cumulative. 


Table  1.  Summary  of  estimated  errors  in  the  budget  retrieval.  Error  percentages  are  valid  for  CK  values  of  about  2.4  X  10  3  and  CK  values 

of  about  1  X  10-3. 


Error  source 

Unmitigated 

error 

Mitigation 

Estimated  error 

Mass  continuity  residual 

>500% 

Use  streamfunction  for  analysis 

0% 

Unresolved  budget 
residual  terms 

±200% 

Constrain  control  volume  size 

±50% 

Shear  flux  term 

-60% 

Estimate  term  using  derived  CD 

±20% 

SST  errors 

—20%  per  1°C  error 

Use  AXBT  data  when  available 

±20% 

Center  errors 

±25%  per  1  km  error  for 
unaveraged  analysis 

Use  high-resolution  radar 

circulation  centers  and  averaging 

±20% 

Analysis  errors 

Large  with  hand  or 
objective  analysis 

SAMURAI 

±20% 

November  2012 


BELL  ET  AL. 


3203 


d.  SAMURAI 

Atmospheric  data  must  be  available  at  regularly 
spaced  intervals  to  numerically  integrate  the  kinematic 
and  thermodynamic  integrals  in  the  budget  calculations. 
For  early  TC  momentum  and  energy  budget  studies  such 
as  that  of  Palmen  and  Riehl  (1957),  data  gridding  was 
done  by  a  manual  map  analysis  that  was  then  inter¬ 
polated  to  regular  intervals.  Manual  analysis  was  largely 
superseded  by  the  development  of  objective  analysis 
techniques  that  created  a  gridded  data  field  by  weighting 
the  observations  by  their  distance  from  a  particular  grid 
point,  such  as  the  Barnes  (1973)  technique  used  in  BM08. 
Modern  data  assimilation  techniques  have  continued  to 
improve  the  analysis  of  data  using  more  sophisticated 
variational  and  probabilistic  techniques.  For  the  cur¬ 
rent  study,  a  two-dimensional  variational  analysis 
technique  called  SAMURAI  was  developed  based 
primarily  on  the  work  of  Ooyama  (1987, 2002)  and  Gao 
et  al.  (2004).  The  SAMURAI  technique  yields  a  maxi¬ 
mum  likelihood  estimate  of  the  atmospheric  state  for 
a  given  set  of  observations  and  error  estimates  by 
minimizing  a  variational  cost  function.  The  technique 
has  several  advantages  over  the  Barnes  (1973)  tech¬ 
nique  used  in  BM08,  including  (i)  observational  error 
specifications  for  different  instrumentation,  (ii)  use  of 
more  complex  observations  such  as  remote  sensing  data, 
(iii)  the  addition  of  balance  constraints  such  as  mass 
continuity,  and  (iv)  a  priori  background  estimates  of  the 
atmospheric  state  when  available. 

A  distinguishing  characteristic  of  the  SAMURAI 
technique  compared  to  other  variational  solvers  is  that 
the  analysis  can  be  performed  directly  in  an  axisymmetric 
cylindrical  coordinate  system.  The  two-dimensional  solver 
improves  the  computational  efficiency  and  minimizes 
potential  errors  in  mass  conservation  that  arise  when 
interpolating  from  a  three-dimensional  domain.  An¬ 
other  distinguishing  characteristic  from  other  varia¬ 
tional  solvers  is  the  use  of  a  Galerkin  approach,  which  is 
similar  to  the  Fourier  spectral  transform  but  uses  the 
cubic  B -spline  as  a  basis  (Ooyama  2002).  The  disad¬ 
vantage  of  the  B -spline  basis  is  that  it  is  not  orthogonal 
and  therefore  requires  an  extra  matrix  to  obtain  the 
spline  coefficients,  but  this  is  a  fair  tradeoff  with  its  other 
desirable  characteristics.  The  basis  is  computationally 
efficient  and  continuously  differentiable  to  second  or¬ 
der,  allowing  for  efficient,  accurate  interpolation  to  ob¬ 
servation  locations,  flexible  incorporation  of  boundary 
conditions,  and  high  numerical  accuracy  of  kinematic 
derivatives.  The  analysis  is  performed  in  a  manner  sim¬ 
ilar  to  the  spectral  transform  method  (Machenhauer 
1979),  transforming  to  and  from  the  spline  coefficients 
and  physical  space  at  each  step  of  the  cost  function 


Date 


Fig.  2.  NHC  best-track  intensity  for  Hurricanes  Fabian  (black 
line)  and  Isabel  (gray  line).  The  six  intensive  observing  periods  used 
in  this  study  are  highlighted  on  2,  3,  4, 12, 13,  and  14  September. 

minimization.  A  more  technical  description  of  SAMURAI 
is  given  in  appendix  B. 

3.  CBLAST  data  analysis 

a.  CBLAST  dataset 

Six  intensive  observing  periods  (IOPs)  were  con¬ 
ducted  in  Hurricane  Fabian  from  2  to  4  September  and 
in  Hurricane  Isabel  from  12  to  14  September  2003  as 
part  of  the  CBLAST  and  National  Oceanic  and  Atmo¬ 
spheric  Administration  (NO AA) /National  Environmen¬ 
tal  Satellite,  Data,  and  Information  Service  (NESDIS) 
Ocean  Winds  experiments  (Black  et  al.  2007).  NHC  best- 
track  intensities  were  estimated  at  greater  than  120  kt 
(62  m  s-1)  during  the  —1600-2300  UTC  time  period  in 
which  observations  were  collected  in  Fabian  and  Isabel 
(Fig.  2),  which  makes  both  storms  major  hurricanes 
during  the  six  IOPs.  Two  NO  A  A  WP-3Ds  (P3s),  the 
NOAA  G-IV,  and  United  States  Air  Force  C-130  aircraft 
collected  in  situ  flight-level  and  dropsonde  observations, 
with  additional  Doppler  radar  and  stepped  frequency 
microwave  radiometer  (SFMR)  data  obtained  by  the  P3s 
only.  The  in  situ  dataset  is  similar  to  that  used  and  de¬ 
scribed  in  detail  in  BM08.  One  of  the  advantages  of 
a  variational  analysis  is  the  ability  to  effectively  combine 
different  observations  based  on  their  individual  error 
characteristics.  The  analysis  composites  were  not  found 
to  be  strongly  sensitive  to  the  magnitude  of  the  prescribed 
errors,  but  the  composites  exhibited  an  unrealistic  level  of 
detail  if  the  errors  were  too  small.  Observation  errors 
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Table  2.  Analysis  times  (UTC)  for  edited  Doppler  radar  data  for  Fabian  during  2-4  September  and  Isabel  12-14  September. 


Fabian 

Isabel 

2  Sep 

3  Sep 

4  Sep 

12  Sep 

13  Sep 

14  Sep 

1712-1726 

1704-1714 

1930-1949 

1650-1659 

1620-1642 

1618-1646 

1930-1941 

1839-1849 

2020-2035 

1721-1726 

1653-1709 

1716-1730 

1953-2002 

1849-1857 

2120-2139 

1828-1846 

1726-1742 

1816-1830 

2014-2027 

1929-1939 

2143-2159 

1901-1910 

1748-1807 

1849-1902 

2101-2114 

1940-1950 

2203-2219 

1943-1948 

1820-1841 

1906-1919 

2115-2123 

2010-2020 

2234-2250 

2006-2019 

1859-1913 

1936-1947 

2206-2219 

2100-2110 

2024-2036 

1916-1930 

1953-2009 

2223-2235 

2114-2124 

2038-2049 

1936-1950 

2013-2029 

2244-2252 

2204-2211 

2115-2130 

1953-2007 

2035-2050 

2026-2039 

2052-2111 

2043-2055 

2112-2125 

2102-2117 

2130-2147 

2151-2207 

2209-2225 

were  set  to  2  m  s-1  for  horizontal  wind  velocities,  2  g  kg-1 
for  water  vapor,  10  g  m-3  for  density,  and  5  kJ  for  moist 
static  energy.  Vertical  velocity  errors  were  set  to  2  m  s_1 
for  flight  level  data  and  4  m  s_1  for  dropsonde  data. 

The  University  of  Massachusetts  deployed  a  SFMR 
for  measuring  surface  wind  speed  and  rain  rate  on  the 
NOAA-42  aircraft.  The  SFMR  is  a  C-band  (~5  cm) 
wavelength,  downward-pointing  radiometer  that  relates 
brightness  temperatures  at  six  different  frequencies  to 
surface  emissivity  and  to  empirically  derived  wind  speeds. 
Since  wind  speed  is  a  nonlinear  observation  operator,  it 
was  assumed  here  that  the  tangential  velocity  comprised 
the  majority  of  the  retrieved  wind  speed.  If  the  surface 
inflow  angle  is  less  than  30°,  then  the  error  caused  by  this 
assumption  is  no  more  than  about  12%.  Given 


uncertainties  in  the  calibration  of  the  instrument  at  very 
high  wind  speeds  and  the  simplification  of  the  obser¬ 
vation  operator,  the  observation  error  was  set  to  a  rel¬ 
atively  high  value  of  10  m  s_1.  Nevertheless,  these 
SFMR  wind  speed  observations  provided  an  important 
constraint  on  the  surface  wind  speed  when  near-surface 
dropsonde  winds  were  not  available. 

Doppler  radar  data  result  in  a  significant  increase  in 
both  azimuthal  and  radial  data  density  for  the  kinematic 
variables  and  add  significant  value  to  the  analysis.  The 
X-band  (~3  cm)  wavelength  tail  Doppler  radar  aboard 
the  NOAA-42  aircraft  employed  a  fore/aft  scanning 
technique  for  all  missions,  which  provided  radial  ve¬ 
locity  data  in  a  cone  about  20°  from  the  track  both  fore 
and  aft  of  the  aircraft.  The  data  were  first  corrected 


Fig.  3.  SST  derived  from  the  TRMM  Microwave  Imager  (average  SST  over  31  August- 
2  September;  color)  and  AXBT  data  released  into  Hurricane  Fabian.  AXBT  numbers  corre¬ 
spond  to  Table  3.  Track  of  Hurricane  Fabian  (dashed  best  track,  2-5  September)  is  shown  for 
reference. 
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Table  3.  SST  observations  (°C)  for  Hurricane  Fabian. 


AXBT  No. 

Time  and  date 

AXBT  SST 

TMI  SST 

SST  difference 

Average  SST 

Analysis  SST 

1 

1725:24  UTC  2  Sep 

26.64 

28.95 

-2.31 

27.795 

2 

2002:29  UTC  2  Sep 

27.61 

28.65 

-1.04 

28.13 

3 

2017:40  UTC  2  Sep 

26.88 

28.5 

-1.62 

27.69 

4 

2049:46  UTC  2  Sep 

26.72 

28.8 

-2.08 

27.76 

27.0 

5 

1712:06  UTC  3  Sep 

27.73 

28.95 

-1.22 

28.34 

6 

1738:25  UTC  3  Sep 

27.4 

28.65 

-1.25 

28.025 

7 

1915:23  UTC  3  Sep 

27.34 

28.65 

-1.31 

27.995 

8 

2034:07  UTC  3  Sep 

27.64 

28.5 

-0.86 

28.07 

27.5 

9 

1944:36  UTC  4  Sep 

27.49 

29.1 

-1.61 

28.295 

28.0 

for  navigational  errors  and  manually  edited  to  remove 
ocean  returns,  radar  side  lobes,  and  other  artifacts  (Oye 
et  al.  1995;  Testud  et  al.  1995;  Bosart  et  al.  2002).  A  large 
number  of  radial  penetrations  were  edited  for  the  cur¬ 
rent  study,  which  included  all  of  the  multiple  dropsonde 
releases.  The  times  of  the  edited  radar  data  used  in 
the  analysis  are  shown  in  Table  2.  Most  of  these  radial 
penetrations  were  used  for  circulation  center  fixes,  but 
a  few  were  too  short  to  obtain  meaningful  center  esti¬ 
mates.  The  Doppler  data  had  an  unambiguous  radial 
velocity  of  less  than  16  m  s-1,  and  frequently  were  aliased 
multiple  times  across  the  Nyquist  interval  because  of  the 
very  high  wind  speeds.  The  Bargen  and  Brown  (1980) 
algorithm  was  used  to  dealias  the  majority  of  the  velocities 
using  the  in  situ  aircraft  flight-level  wind  as  a  reference, 
but  gaps  and  noise  in  the  data  required  additional 
manual  unfolding  of  many  radar  beams.  The  Doppler 


Fig.  4.  Storm-relative  AXBT  splash  locations.  AXBT  numbers 
correspond  to  Fig.  3.  Gray  annulus  corresponds  to  approximate 
eyewall  location  at  30-km  radius  from  the  storm  center. 


velocities  were  then  averaged  along  each  beam  with 
a  spatial  resolution  that  approximately  matched  the  di¬ 
ameter  of  the  beam  with  increasing  range,  in  order  to 
reduce  noise  and  the  voluminous  number  of  individual 
radar  observations.  The  Doppler  velocity  error  was  as¬ 
signed  as  the  sum  of  the  spectrum  width  of  the  velocity 
measurements  (Keeler  and  Ellis  2000)  plus  an  estimated 
error  in  the  terminal  fall  speed  of  the  precipitation,  with 
a  minimum  error  of  2  m  s-1. 

b.  Sea  surface  temperature 

Direct  measurements  of  the  SST  were  made  during 
the  Hurricane  Fabian  missions  by  AXBTs  released  by 
the  NO  A  A  aircraft.  A  comparison  of  the  AXBT  tem¬ 
peratures  and  Tropical  Rainfall  Measuring  Mission 
(TRMM)  Microwave  Imager  (TMI)  satellite  estimates 
from  the  3-day  period  prior  to  the  first  mission  is  illus¬ 
trated  in  Fig.  3  and  summarized  in  Table  3.  The  TMI- 
derived  sea  surface  temperatures  are  near  28.5°C  along 
the  track  throughout  the  analysis  domain,  but  the  AXBT 
temperatures  are  consistently  1°-2°C  lower  at  the  same 
locations.  Since  decreases  of  the  SST  are  known  to  occur 
during  and  after  the  storm  passage,  the  splash  locations 
of  the  measurements  were  plotted  in  a  storm-relative 


Fig.  5.  NHC  best  tracks  for  Hurricanes  Fabian  (black  line)  and 
Isabel  (gray  line)  and  radar-derived  tracks  (color  lines  for  Hurri¬ 
canes  Fabian  and  Isabel).  Color  lines  indicate  the  analysis  period 
for  each  aircraft  mission  used  in  this  study  (see  inset). 
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Fig.  6.  Spatial  distribution  of  observations  from  2  September  Hurricane  Fabian  missions.  (a),(c)  Nonradar  and  (b), 
(d)  radar  observations  are  plotted  in  the  (a),(b)  radius-height  and  (c),(d)  radius-azimuth  planes.  Red  triangles 
indicate  dropsonde  observations,  blue  circles  indicate  flight-level  in  situ  observations,  green  squares  indicate  SFMR 
observations,  and  black  dots  indicate  Doppler  radar  observations. 


coordinate  system  (Fig.  4).  All  AXBT  measurements 
were  near  the  eyewall  (gray  annulus)  with  the  exception 
of  AXBTs  4  and  8.  On  2  September,  the  two  AXBTs  on 
the  right  side  of  the  eyewall  are  about  1°C  lower  than  the 
one  released  on  the  left  side.  Given  the  relatively  good 
azimuthal  sampling  around  the  eyewall,  an  average  value 
of  27°C  was  chosen  for  this  day.  On  3  September,  all  of 
the  AXBTs  near  the  eyewall  agree  to  within  0.4°C  and 
are  consistently  1.2°C  lower  than  the  TMI  temperatures. 
A  slightly  higher  value  of  27.5°C  was  utilized  on  this  day. 
It  is  interesting  that  the  general  location  of  Fabian  on  3 
September  coincides  with  that  of  Isabel  on  13  September. 
Since  the  SST  is  estimated  at  27.5°C  on  both  of  these  days, 
it  appears  that  some  of  the  ocean  cooling  that  occurred 
during  Fabian’s  passage  recovered  in  the  10  days  be¬ 
tween  the  two  storms.  Only  a  single  AXBT  was  available 
on  4  September,  which  was  adjusted  upward  slightly  to 
28°C  given  the  warmer  SSTs  derived  from  TMI. 


A  comprehensive  analysis  of  the  SST  for  the  potential 
intensity  estimates  for  Hurricane  Isabel  was  performed 
by  BM08  using  satellite,  buoy,  and  airborne  radiometer 
measurements.  Comparisons  of  the  3-day  TMI  micro- 
wave-derived  product  with  the  Advanced  Very  High 
Resolution  Radiometer  (AVHRR)  infrared  estimates 
are  good,  but  given  the  high  bias  of  the  TMI  estimates  in 
the  Fabian  case,  the  29°C  on  14  September  used  in  BM08 
was  adjusted  downward  slightly  here  to  28.5°C.  The 
28.5°  and  27.5°C  sea  surface  temperatures  obtained 
previously  for  12  and  13  September,  respectively,  were 
not  modified. 

c.  Circulation  centers 

There  are  a  variety  of  ways  to  define  the  center  of 
a  tropical  cyclone,  but  an  optimal  center  for  the  current 
study  is  one  that  minimizes  aliasing  associated  with  a  trun¬ 
cated  Fourier  representation  of  the  angular  momentum 
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Fig.  7.  Spatial  distribution  of  nonradar  observations  from  Hurricane  Fabian  missions  on  (a),(b)  3  and  (c),(d)  4  September. 
Observations  are  plotted  in  the  (a),(c)  radius-height  and  (b),(d)  radius-azimuth  planes.  Symbols  are  as  in  Fig.  6. 


and  total  energy  in  azimuthal  wavenumbers.  The  “cir¬ 
culation  center”  proposed  by  Marks  et  al.  (1992)  meets 
this  criterion  and  is  defined  as  the  center  that  maximizes 
the  axisymmetric  tangential  wind  at  the  radius  of  maxi¬ 
mum  wind  (RMW)  derived  from  Doppler  radar  data. 
Since  there  is  an  inflection  point  in  the  radial  gradient  of 
tangential  wind  at  the  RMW,  the  circulation  center  also 
maximizes  the  vorticity  at  the  RMW  and  the  circulation 
inside  the  RMW.  The  circulation  center  maximizes  the 
symmetric  tangential  wind  by  definition,  and  therefore 
minimizes  aliasing  onto  higher-order  azimuthal  wave- 
numbers.  The  use  of  a  broad  swath  of  Doppler  radar- 
derived  wind  data  also  helps  to  eliminate  mesovortex 
centers  that  could  be  used  mistakenly  by  analyzing  only 
in  situ  data  (Willoughby  and  Chelmow  1982). 

To  determine  the  circulation  center,  a  two-dimensional 
“simplex”  search  (Nelder  and  Mead  1965)  was  used  on  a 


1.5-km  gridded  dual-Doppler  horizontal  wind  field 
(Mohr  et  al.  1986)  to  find  the  point  that  maximized 
the  azimuthally  averaged  tangential  wind  in  an  an¬ 
nulus  at  a  specified  radius.  Centers  were  deduced  at 
1-km  altitude  for  each  radial  penetration  of  the  air¬ 
borne  radar  and  interpolated  over  time  using  a  cubic 
spline  to  create  a  high-resolution  track.  The  track  was 
reviewed  and  adjusted  subjectively  to  remove  out¬ 
liers  and  centers  biased  by  local  wind  maxima,  with 
the  accuracy  of  the  derived  high-resolution  track  es¬ 
timated  at  about  3  km.  The  NHC  best  track  and  the 
high-resolution  tracks  derived  from  the  simplex 
searches  are  shown  for  all  six  missions  in  Fig.  5.  As 
expected,  the  tracks  are  nearly  coincident,  but  the 
higher-resolution  motion  is  apparent  in  the  colored  tracks. 
The  largest  discrepancy  appears  to  be  on  13  September, 
where  the  differences  are  likely  due  to  the  different 
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Fig.  8.  Spatial  distribution  of  nonradar  observations  from  Hurricane  Isabel  on  (a),(b)  12,  (c),(d)  13,  and  (e),(f)  14 
September.  Observations  are  plotted  in  the  (a),(c),(e)  radius-height  and  (b)(d),(f)  radius-azimuth  planes.  Symbols 
are  as  in  Fig.  6. 
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Fig.  9.  Hurricane  Fabian  axisymmetric  tangential  wind  (color),  radial  wind  (5  m  s-1  con¬ 
tours),  and  secondary  circulation  (vectors)  from  2  to  4  September.  Dashed  contours  indicate 
inflow  and  solid  contours  indicate  outflow. 


center  definitions,  but  are  still  well  within  the  large  eye 
on  this  day. 

d.  SAMURAI 

The  SAMURAI  technique  was  conducted  for  the 
six  CBLAST  missions  using  the  data  sources  and  cir¬ 
culation  centers  described  above.  The  radar  and  non¬ 
radar  observation  data  distribution  for  the  2  September 
mission  is  shown  in  Fig.  6  in  both  the  radius-height  and 
polar  planes.  It  is  clear  that  the  radar  data  dominate  the 
spatial  distribution  because  of  the  wide  swath  of  the  tail 
Doppler  radar.  Comprehensive  radar  coverage  provides 
good  confidence  in  the  kinematic  fields  throughout  the 
domain,  except  near  the  eye  where  scatterers  are  limited. 
Excellent  dropsonde  coverage  is  also  obtained  in  the  main 
eyewall  region  near  the  RMW,  with  bracketing  observa¬ 
tions  at  flight  level  and  at  the  surface  from  the  SFMR. 
Good  azimuthal  sampling  was  achieved  by  the  multiple 
sequences  across  the  eyewall.  Since  only  the  flight  level 


and  dropsonde  data  provide  thermodynamic  data,  the 
highest-quality  energy  analysis  is  limited  to  the  —20-50- 
km  annulus.  Since  there  are  no  high  confidence  a  priori 
estimates  of  the  structure  available,  other  than  the  pre¬ 
vious  objective  analysis  with  the  same  dataset,  the  data  is 
only  weakly  constrained  in  poor  data  regions.  By  design 
of  the  CBLAST  experiment,  the  data-rich  regions  co¬ 
incide  with  the  regions  of  highest  wind  speeds  that  are  of 
interest  for  the  current  study. 

The  data  distributions  of  nonradar  observations  for 
the  remainder  of  the  missions  are  shown  in  Figs.  7  and  8. 
Radar  sampling  was  similar  on  all  days  to  the  first  mis¬ 
sion  and  is  not  shown.  The  dropsonde  distributions  are 
similar  to  the  first  mission,  with  limited  data  outside  of 
40-km  radius  but  excellent  azimuthal  sampling  on  3 
September,  and  limited  azimuthal  sampling  on  4  Sep¬ 
tember.  All  of  the  Isabel  missions  are  characterized  by 
comprehensive  radial  and  azimuthal  sampling  (Fig.  8). 
Note  that  distributions  are  slightly  different  from 
those  shown  in  BM08  due  to  the  use  of  radar-derived 
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Fig.  10.  Hurricane  Fabian  axisymmetric  angular  momentum  (color),  total  energy  (contours), 
and  secondary  circulation  (vectors)  from  2  to  4  September.  The  thick  contour  is  350  kJ  kg-1;  solid 
lines  and  dashed  lines  indicate  values  above  and  below  350  kJ  kg-1,  respectively,  with  a  contour 
interval  of  1  kJ  kg-1. 


circulation  centers  instead  of  flight-level  winds  as  de¬ 
scribed  in  section  3b,  but  the  patterns  are  qualitatively 
the  same.  The  number  of  observations  included  in  each 
composite  varies  from  about  500  000  to  more  than  1.7 
million  on  14  September  (due  to  the  large  domain  re¬ 
quired  to  encompass  the  45-km  RMW).  The  difference 
between  the  observations  and  analysis  yielded  an  ap¬ 
proximately  normal  distribution  with  a  mean  and  median 
near  zero,  linear  correlations  between  the  observations 
and  analysis  near  0.99,  and  linear  slope  and  bias  values 
near  one  and  zero,  respectively,  for  all  days  (not  shown). 
The  statistics  indicate  that  the  SAMURAI  technique 
had  high  fidelity  to  the  observations  and  suggest  that  the 
TCs  were  very  azimuthally  symmetric  during  all  six 
missions. 

The  derived  wind  fields  for  Hurricane  Fabian  are 
shown  in  Fig.  9.  The  tangential  winds  depict  a  general 
weakening  over  the  three  analysis  days,  with  a  decrease 
in  the  depth  of  the  strongest  winds  and  slow  weakening 


of  the  maximum  tangential  wind  by  the  third  day.  The 
reduction  in  the  depth  and  intensity  of  the  tangential 
wind  was  concurrent  with  a  reduction  in  the  inflow  depth 
and  magnitude  over  the  3  days.  The  primary  updraft  ap¬ 
pears  to  have  weakened  also  during  these  3  days.  The 
RMW  remained  consistently  just  inside  the  30-km  radius, 
with  a  sharp  gradient  of  tangential  wind  toward  the 
center.  Note  that  the  low-level  inflow  continues  past  the 
RMW  on  all  3  days,  with  the  flow  turning  upward  near 
the  high  tangential  wind  gradient  region.  The  absolute 
angular  momentum  M  and  total  energy  E  are  shown  in 
Fig.  10.  The  secondary  circulation  generally  follows  the 
contours  of  M  above  the  inflow  layer,  with  more  radi¬ 
ally  tilted  M  surfaces  outside  the  RMW  on  the  second 
and  third  days.  The  E  contours  depict  also  a  toroidal 
shape  that  is  similar  to,  but  not  congruent  with,  the  M 
contours  and  secondary  circulation.  Note  that  the  closed 
E  contours  at  the  outer  radii  (~50-km  radius)  are  likely 
artifacts  of  the  lack  of  data  where  the  analysis  relaxes 
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Fig.  11.  Hurricane  Isabel  axisymmetric  tangential  wind  (color),  radial  wind  (5  m  s-1  contours), 
and  secondary  circulation  (vectors)  from  12  to  14  September. 


back  to  the  prescribed  reference  state.  Similarly,  the 
thermodynamic  structure  inside  of  20-km  radius  may  be 
unreliable  because  of  a  data  void.  A  reduced  radial 
energy  gradient  is  apparent  on  each  consecutive  day, 
which  is  consistent  with  the  weakening  trend  in  the 
kinematic  variables. 

The  analyzed  kinematic  fields  for  Hurricane  Isabel 
for  all  three  missions  are  shown  in  Fig.  11.  Note  that  the 
inner  radius  shown  is  15  km  for  clarity  because  of  the 
larger  domain,  although  the  analysis  domain  was  ex¬ 
tended  to  the  TC  center.  The  qualitative  features  of  the 
SAMURAI  technique  are  similar  to  the  Barnes  analysis 
used  in  BM08,  but  with  enforced  mass  continuity  of  the 
secondary  circulation  and  a  notably  sharper  resolution 
of  the  boundary  layer  inflow.  The  characteristic  “over¬ 
shoot”  of  the  inflow  past  the  RMW  is  apparent  during 
all  three  missions  as  was  the  case  in  the  Fabian  analysis, 
with  a  sloping  inflow  top  that  coincides  with  the  height 
of  the  maximum  tangential  wind.  A  gradual  weakening 
and  expansion  of  the  tangential  wind  as  described  in 
BM08  is  evident,  with  a  well-defined  outflow  and  up¬ 
draft  core  just  above  the  maximum  tangential  wind  on 
all  3  days.  The  absolute  angular  momentum  and  total 
energy  analyses  for  Hurricane  Isabel  are  shown  in 
Fig.  12.  The  consistent  structures  of  the  secondary  cir¬ 
culation,  M,  and  E  contours  are  similar  to  those  in  the 


analyses  of  Fabian.  The  total  energy  structure  is  also 
qualitatively  similar  to  the  moist  entropy  structure  shown 
by  BM08  (cf.  their  Fig.  5),  as  would  be  expected  for 
similar  thermodynamic  variables. 

The  kinematic  and  thermodynamic  structures  appar¬ 
ent  in  the  six  analyses  are  consistent  in  their  depiction  of 
gradually  weakening  storms  just  past  their  peak  intensity. 
The  two  different  TCs  and  six  different  missions  provide 
a  good  dataset  for  the  flux  retrievals,  since  the  TCs  were 
both  in  an  approximate  steady  state  but  span  20  m  s”1 
in  their  peak  tangential  wind  speeds  and  have  slightly 
different  RMW,  total  energy,  and  secondary  circulation 
characteristics. 

4.  Air-sea  fluxes 

The  absolute  angular  momentum  and  total  energy 
budgets  were  resolved  using  the  SAMURAI-derived 
axisymmetric  composites  shown  in  the  previous  section. 
A  variety  of  control  volumes  was  used  to  take  into  ac¬ 
count  the  errors  associated  with  the  unresolved  budget 
terms  and  quantitatively  estimate  the  uncertainty  in  the 
derived  fluxes.  To  determine  the  top  of  the  control 
volume,  traditional  formalism  for  the  planetary  bound¬ 
ary  layer  assumes  that  the  magnitude  of  the  turbulent 
fluxes  decreases  to  zero  at  the  top  of  the  mixed  layer, 
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Fig.  12.  Hurricane  Isabel  axisymmetric  angular  momentum  (color),  total  energy  (contours), 
and  secondary  circulation  (vectors)  from  12  to  14  September.  The  thick  contour  is  350  kJ  kg-1; 
solid  lines  and  dashed  lines  indicate  values  above  and  below  350  kJ  kg-1,  respectively,  with 
a  contour  interval  of  1  kJ  kg-1. 


but  the  TCBL  has  a  distinct  separation  of  the  boundary 
layer  height  defined  mechanically  and  thermodynami¬ 
cally  (Montgomery  et  al.  2006;  Smith  and  Montgomery 
2009;  Zhang  et  al.  2011).  In  between  the  outer  rain- 
bands  of  Hurricanes  Fabian  and  Isabel,  the  vertical 
momentum  and  humidity  fluxes  decreased  to  zero  at 


700  m,  but  the  virtual  potential  temperature  mixed 
layer  depth  was  400  m,  and  the  inflow  layer  extended  to 
1  km  (Zhang  et  al.  2009).  Mixed  layer  depths  and  in¬ 
flow  heights  are  similar  in  the  eyewall  (Montgomery 
et  al.  2006),  but  the  TCBL  may  have  different  turbu¬ 
lence  heights  in  the  main  updraft  region  compared  with 
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Fig.  13.  Control  volumes  used  for  flux  retrievals  on  12  September.  Tangential  wind  (color)  and 
secondary  circulation  (vectors)  are  shown  for  reference.  Thick  white  lines  indicate  the  borders  of 
overlapping  control  volumes.  Small  apparent  boxes  near  the  corners  of  the  control  volumes  are 
artifacts  resulting  from  the  overlap.  An  example  control  volume  is  highlighted  in  black. 
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the  inflow  region  outside  the  eyewall  (Smith  et  al.  2008). 
Sensitivity  tests  suggested  that  volume  tops  below  and 
near  the  top  of  the  inflow  layer  were  the  most  appropriate 
(appendix  A),  and  a  range  of  control  volume  depths 
from  400  to  1000  m  was  therefore  utilized.  Results  with 
a  fixed  1-km  control  volume  top  were  within  10%  and 
20%  of  the  mean  CD  and  CK  values  over  all  volume 
depths,  respectively. 

The  bottom  of  the  control  volume  was  set  within  the 
surface  layer  at  either  10-  or  100-m  altitude.  Variable 
control  volume  widths  from  10  to  22  km  were  prescribed 
in  2-km  increments,  and  the  inner  radius  of  the  volume 
was  varied  by  2-km  increments  around  a  central  value  of 
80%  of  the  RMW.  These  volumes  spanned  the  region 
in  and  around  the  eyewall,  such  that  the  average  wind 
speed  in  the  volume  was  always  above  50  m  s”1  for  every 
mission.  The  aspect  ratio  (width  divided  by  height)  of  all 
control  volumes  was  between  20  and  30  as  suggested  by 
the  sensitivity  analysis.  An  illustration  of  the  control  vol¬ 
umes  used  on  12  September  is  shown  in  Fig.  13.  The  72 
white  control  volumes  overlap  throughout  the  eyewall 
region,  with  one  of  the  example  control  volumes  high¬ 
lighted  in  black  for  clarity.  Using  multiple  control  vol¬ 
umes  yielded  a  total  of  72  flux  “samples”  per  mission, 
representing  a  range  of  volume  shapes,  sizes,  and  loca¬ 
tions  centered  around  the  RMW.  The  spectrum  of  control 
volumes  adequately  represents  a  reasonable  statistical 
sample  for  obtaining  a  mean  and  standard  deviation  of  the 
budget  residual  at  the  eyewall  on  each  day. 

The  retrieved  surface  stress  values  are  shown  in  Fig. 
14a.  Although  considerable  spread  exists  in  the  individual 
samples,  a  general  agreement  in  the  stress  magnitudes  is 
found  with  a  mean  value  of  9.4  N  m”2  and  standard  de¬ 
viation  of  4.6  N  m-2.  A  linear  fit  of  the  data  has  a  corre¬ 
lation  coefficient  of  0.22  and  only  a  slight  upward  trend 
with  wind  speed.  Missions  on  2  and  13  September  have 
the  highest  stress  but  also  two  of  the  deepest  inflow  layers. 
The  14  September  analysis  has  a  deep  inflow  layer  also, 
but  with  a  stress  more  similar  to  the  other  missions.  The 
retrieved  stress  on  12  September  is  low  compared  to  the 
other  missions.  It  is  not  clear  that  the  low  stress  is  due  to 
physical  differences  on  this  day;  it  seems  more  likely  that 
the  scatter  is  due  to  inherent  uncertainties  in  the  meth¬ 
odology.  Removing  this  mission  from  the  sample  im¬ 
proves  the  correlation  coefficient  to  0.53  and  results  in 
a  more  distinct  upward  trend  with  wind  speed.  The  de¬ 
rived  friction  velocity  is  shown  in  Fig.  14b.  The  scatter  is 
less  than  in  the  stress  retrieval  due  to  the  square  root 
dependence  on  the  stress,  with  a  mean  value  of  2.8  ms-1 
and  a  0.7  m  s”1  standard  deviation.  These  estimates  in¬ 
dicate  also  a  slight  upward  trend  with  wind  speed  with 
a  correlation  of  0.24,  or  0.53  with  the  12  September  mis¬ 
sion  removed. 


Fig.  14.  (a)  Derived  surface  stress  and  (b)  derived  friction  velocity 
from  the  absolute  angular  momentum  budget  retrievals.  Gray  dots 
indicate  individual  samples  from  different  control  volumes,  and 
large  symbols  indicate  mean  values  from  each  research  mission. 
Error  bars  indicate  one  standard  deviation  in  derived  quantity 
(vertical)  and  average  surface  wind  speed  (horizontal). 

The  retrieved  surface  enthalpy  flux  is  shown  in  Fig.  15. 
These  are  the  first  known  estimates  of  the  enthalpy  flux 
at  wind  speeds  greater  than  50  m  s-1,  and  likely  the 
highest  estimates  of  heat  flux  recorded  over  the  ocean 
surface.  There  is  less  scatter  than  in  the  stress  retrieval, 
and  an  increasing  linear  correlation  with  wind  speed  of 
0.81.  The  increasing  trend  with  wind  speed  is  larger  in 
the  individual  retrievals  for  each  mission  than  in  the 
mean,  which  may  be  due  to  the  use  of  thinner  control 
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and  total  energy  budget  retrievals.  Symbols  are  as  in  Fig.  14.  (b)  Enthalpy  Exchange  Coefficient  (CK) 


volumes  for  the  higher  wind  speeds.  The  wind  speed 
dependence  of  the  mean  enthalpy  flux  from  each  mission 
increases  from  764  W  m~2  at  52  m  s”1  to  2189  W  m-2  at 
72  m  s_1.  The  reasons  for  the  reduced  scatter  and  better 
correlation  in  the  energy  budget  compared  to  the  mo¬ 
mentum  budget  are  not  apparent.  It  would  appear  that 
the  derived  thermodynamic  structure  was  steadier  on  the 
~6-h-mission  time  scale  than  the  kinematic  structure,  but  it 
is  impossible  to  estimate  the  thermodynamic  tendency 
with  the  current  dataset. 

The  derived  10-m  drag  and  enthalpy  exchange  co¬ 
efficients  are  shown  in  Fig.  16.  The  CK  estimates  show 
better  agreement  than  the  CD  estimates,  as  would  be 
expected  from  better  agreement  of  the  enthalpy  flux 
retrievals.  The  mean  CD  is  2.4  X  10-3  with  a  standard 
deviation  of  1.1  X  10-3,  which  corresponds  to  about 
46%  uncertainty  at  the  67%  confidence  level  and  about 
93%  uncertainty  at  the  95%  confidence  level.  The  mean 
CK  is  1.0  X  10-3  with  a  standard  deviation  of  0.4  X  10-3, 
which  corresponds  to  about  40%  uncertainty  at  the  67% 
confidence  level  and  about  80%  uncertainty  at  the  95% 
confidence  level.  These  percentage  uncertainties  are  con¬ 
sistent  with  a  sensitivity  analysis  described  in  appendix  A 
that  indicated  an  approximately  50%  standard  deviation 
associated  with  the  neglected  budget  terms.  The  CD  results 
do  not  indicate  a  significant  dependence  on  wind  speed, 
with  a  slightly  decreasing  linear  correlation  of  only  0.11. 
The  CK  results  show  a  slight  increase  with  wind  speed,  with 
a  slope  of  0.03  (m  s-1)-1  and  a  higher  linear  correlation  of 
0.53.  Given  the  uncertainties  and  the  weak  slopes  of  both 
fits,  there  is  no  statistically  significant  change  in  either  CD 
or  CK  in  the  wind  speed  range  analyzed. 


Fig.  16.  Derived  10-m  bulk  exchange  coefficients  from  the  abso¬ 
lute  angular  momentum  and  total  energy  budget  retrievals.  Symbols 
are  as  in  Fig.  14. 


The  histograms  of  the  CD  and  CK  samples  shown  in 
Fig.  17  suggest  unimodal  statistical  distributions  with 
some  spread  attributable  to  random  errors.  The  statis¬ 
tical  distributions  are  slightly  skewed  to  the  right  (0.56 
and  0.45  skewness  for  CD  and  CK,  respectively),  sug¬ 
gesting  that  the  mean  has  a  slight  low  bias  compared  to  a 
normal  distribution.  Negative  kurtosis  (-0.58  and  -0.12, 
respectively)  indicates  a  rounder  peak  and  shorter  tails 
than  a  normal  distribution.  The  cumulative  probability 
distribution  functions  (CDFs;  Fig.  17b)  graphically 
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Fig.  17.  (a)  Histogram  of  the  percentage  of  samples  in  0.5  X  10-3 
bins  for  CD  (black)  and  CK  (gray)  for  all  retrievals,  (b)  Cumulative 
probability  distributions  for  CD  (black)  and  CK  (gray)  retrievals. 


illustrate  the  probability  that  the  exchange  coefficients 
are  at  or  below  a  particular  magnitude.  These  probabili¬ 
ties  assume  that  the  432  samples  obtained  from  varying 
the  control  volume  over  each  of  the  six  missions  accu¬ 
rately  represent  samples  from  the  true  population  dis¬ 
tribution  of  exchange  coefficient  magnitudes.  This 
assumption  is  based  on  the  expectation  that  the  varia¬ 
tion  of  control  volumes  reasonably  depicts  the  un¬ 
certainty  associated  with  the  unresolved  budget  terms, 
and  that  the  remaining  errors  are  random.  The  CDF  for 


(a)  Derived  CK/  CD  ratio 
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Fig.  18.  (a)  Ratio  of  CKICD  from  budget  retrievals.  Symbols  are  as  in 
Fig.  14.  (b)  Cumulative  probability  distribution  for  CKICD  retrievals. 


CD  has  a  steep  slope  above  1  X  10-3,  with  50%,  70%,  and 
90%  probabilities  that  CD  is  less  than  2,  3,  and  4  X  10-3, 
respectively.  The  maximum  estimated  CD  is  less  than  5  X 
10-3  with  99%  probability,  but  this  upper  limit  is  not  that 
useful  given  prior  estimates  and  physical  expectations. 
The  slope  of  the  CDF  for  CK  is  much  flatter  than  that  of 
the  corresponding  slope  for  CD.  The  maximum  estimated 
CK  is  less  than  2  X  10-3  with  99%  probability,  with  a  60% 
probability  that  the  magnitude  is  less  than  1  X  10 _3. 

The  ratios  of  CkICD  calculated  from  the  different 
samples  and  the  cumulative  probability  distribution  are 
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Fig.  19.  Wind  speed  dependence  of  CD  from  this  study  (green  circles)  compared  with  pre¬ 
vious  studies.  Black  symbols  indicate  data  adapted  from  French  et  al.  (2007)  and  blue  symbols 
indicate  data  adapted  from  Vickery  et  al.  (2009).  Red  line  indicates  measured  (thick)  and 
extrapolated  (thin)  Large  and  Pond  (1981)  drag  coefficient. 


shown  in  Fig.  18.  With  the  exception  of  12  September, 
the  retrieved  ratios  are  in  relatively  good  agreement. 
Dividing  the  mean  CkICD  gives  an  average  ratio  of  0.4, 
with  a  variation  from  0.17  to  1.05  by  adding  and  sub¬ 
tracting  a  single  standard  deviation  of  each  coefficient  in 
each  direction.  The  mean  value  of  all  ratios  calculated 
individually  from  each  sample  is  0.48  with  a  standard 
deviation  of  0.27.  A  linear  fit  shows  an  increasing  ratio 
with  wind  speed  with  a  correlation  of  0.58,  but  the  in¬ 
crease  is  primarily  due  to  the  large  ratios  above  70  m  s-1. 
It  is  unclear  why  the  12  September  mission  appears  to 
be  an  outlier,  but  this  is  probably  due  to  unresolved 
budget  terms  as  opposed  to  a  significant  change  in  the 
ratio  above  70  m  s-1.  It  cannot  be  ruled  out  that  CD  is 
reduced  at  these  extreme  wind  speeds,  but  since  Ck  on 
12  September  agrees  well  with  the  other  five  cases,  CD 
derived  for  this  day  seems  low.  This  discrepancy  could 
be  due  to  a  potential  eyewall  replacement  shortly  after 
the  aircraft  mission  (BM08),  indicating  a  larger  un¬ 
resolved  tendency  term  than  on  the  other  days.  Even 
with  the  outlier  ratio  included,  the  cumulative  proba¬ 
bility  distribution  of  all  CKICD  ratios  shown  in  Fig.  19b 
indicates  an  80%  probability  that  the  ratio  is  less  than 
0.75,  and  a  93%  probability  that  it  is  less  than  1.0. 
Without  the  12  September  ratio  these  probabilities  in¬ 
crease  to  91%  and  97%,  respectively  (not  shown). 

5.  Summary  and  discussion 

A  budget  methodology  to  retrieve  air-sea  fluxes  of 
momentum  and  enthalpy  was  performed  using  data 
from  six  CBLAST  missions  into  major  hurricanes  Fabian 
and  Isabel  (2003)  at  a  range  of  surface  wind  speeds  from 
52  to  72  m  s_1.  This  study  presents  what  are  believed  to 
be  the  first  estimates  of  enthalpy  fluxes  ( CK )  and  the 
ratio  of  CK/CD  in  major  hurricanes  at  wind  speeds 
greater  than  50  m  s”1.  The  results  presented  here  also 
augment  the  limited  surface  stress  and  CD  estimates  at 
these  wind  speeds.  Several  potential  errors  were 


identified  by  sensitivity  analyses  using  simulated  data 
and  the  errors  were  optimally  mitigated  to  reduce  the 
uncertainty  of  the  estimates.  The  main  sources  of  un¬ 
certainty  identified  were  (i)  unresolved  budget  terms,  (ii) 
sea  surface  temperature  errors,  (iii)  circulation  center 
errors,  and  (iv)  analysis  errors.  The  “shear  flux”  term, 
representing  the  interaction  of  turbulent  momentum 
fluxes  and  the  mean  shearing  flow,  was  found  to  be  im¬ 
portant  in  the  total  energy  budget  and  was  included 
through  an  estimate  of  CD  obtained  from  the  absolute 
angular  momentum  budget.  Quantitative  uncertainty 
estimates  derived  by  varying  the  budget  control  volumes 
within  the  six  different  analyses  are  consistent  with  un¬ 
certainties  estimated  from  sensitivity  tests  using  simulated 
observations. 

The  mean  CD  estimate  for  wind  speeds  above 
50  m  s”1  from  this  study  is  2.4  X  10-3  with  a  standard 
deviation  of  1.1  X  10-3.  Estimates  of  CD  from  the  cur¬ 
rent  research  are  shown  in  Fig.  19  in  relation  to  previous 
studies.  The  black  symbols  are  adapted  from  French 
et  al.  (2007),  where  the  black  circles  along  the  thick 
black  line  indicate  the  CBLAST  estimates  from  that 
study,  squares  along  the  thin  line  indicate  estimates  from 
Powell  et  al.  (2003),  and  diamonds  along  the  thin  dashed 
line  indicate  laboratory  estimates  from  Donelan  et  al. 
(2004).  The  blue  symbols  are  adapted  from  Vickery 
et  al.  (2009)  showing  the  flux-profile  estimates  obtained 
from  dropsondes.  The  current  results  are  shown  in 
green,  along  with  95%  confidence  intervals  in  black.  The 
CD  estimates  from  this  study  are  in  general  agreement 
with  the  previous  studies  that  CD  does  not  continue  to 
increase  beyond  about  30  m  s”1,  though  some  increase 
cannot  be  strictly  ruled  out  given  the  uncertainties  in  the 
estimates.  At  the  95%  confidence  level,  the  uncertainty 
does  not  preclude  the  possibility  that  the  drag  coefficient 
is  near  the  extrapolated  Large  and  Pond  (1981)  curve 
(dash-dotted  curve),  but  the  cumulative  probability  dis¬ 
tribution  indicates  that  the  magnitude  has  a  90%  prob¬ 
ability  of  being  less  than  4  X  10-3. 
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Fig.  20.  Wind  speed  dependence  of  CK  from  this  study  (green  squares)  compared  with 
previous  studies.  ASIST  laboratory  results  (blue  circles)  and  CBLAST  (red  triangles)  mea¬ 
surements  shown  with  HEXOS  results  (gray  crosses)  adapted  from  Haus  et  al.  (2010).  The 
mean  and  95%  confidence  intervals  are  shown  in  black. 


The  mean  CK  estimate  for  wind  speeds  above  50  m  s” 
from  this  study  is  1.0  X  10-3  with  a  standard  deviation  of 
0.4  X  10“ 3.  Estimates  of  CK  from  the  current  research 
are  shown  in  Fig.  20  with  previous  estimates.  The  cur¬ 
rent  results  are  shown  in  green,  along  with  Humidity 
Exchange  over  the  Sea  (HEXOS;  DeCosmo  et  al.  1996; 
gray  crosses),  CBLAST  (Zhang  et  al.  2008;  red  trian¬ 
gles),  and  laboratory  results  (blue  circles)  adapted  from 
Haus  et  al.  (2010).  The  mean  and  95%  confidence  in¬ 
tervals  for  the  HEXOS,  CBLAST,  and  current  results 
are  indicated  with  black  lines.  The  uncertainty  at  the 
95%  confidence  level  is  higher  in  the  current  study  com¬ 
pared  to  results  from  the  other  field  experiments,  but  the 
envelope  of  CK  magnitudes  is  similar  in  all  the  field  es¬ 
timates.  The  uncertainties  in  CK  and  CD  are  similar  in 
terms  of  a  percentage  error,  but  CK  has  a  lower  un¬ 
certainty  in  terms  of  absolute  error.  There  is  a  good 
linear  correlation  of  the  enthalpy  flux  with  wind  speed 
(R  =  0.81;  see  Fig.  15)  and  the  derived  CK  is  consistent 
across  the  different  aircraft  missions.  A  linear  fit  of  CK 
indicates  a  slight  increase  with  wind  speed  above  50  m  s”1, 
but  the  increase  is  well  within  the  uncertainty  range  and 


cannot  be  concluded  definitively.  The  cumulative  prob¬ 
ability  distribution  indicates  that  CK  has  an  approxi¬ 
mately  70%  probability  of  being  equal  to  or  less  than 
1.2  X  10-3,  which  is  the  approximate  value  determined  at 
wind  speeds  greater  than  15  m  s_1.  These  results  suggest 
that  it  is  probable  that  the  magnitude  of  CK  is  not  de¬ 
pendent  on  wind  speed  in  major  hurricane  conditions. 
Since  the  current  results  implicitly  include  the  effects  of 
sea  spray,  these  results  would  also  suggest  that  spray  ef¬ 
fects  do  not  change  CK  significantly.  One  hypothesis  is 
that  the  spray  flux  has  simply  replaced  the  interfacial 
flux  at  these  wind  speeds  (Haus  et  al.  2010),  but  the  net 
enthalpy  flux  is  similar.  These  results  suggest  also  that  sea 
spray  parameterizations  that  increase  CK  at  high  wind 
speeds  may  need  revision  (Bao  et  al.  2011;  Andreas  2011). 
However,  a  partition  between  the  spray  and  interfacial 
fluxes  cannot  be  assessed  by  the  current  methodology. 

The  retrieved  CK/CD  ratios  from  the  current  research 
are  shown  with  previous  estimates  in  Fig.  21.  These  re¬ 
sults  indicate  that  the  ratio  is  likely  less  than  the  1.0 
estimate  derived  by  Emanuel  (2003)  at  the  ~28°C 
temperature  range,  and  perhaps  may  be  as  low  as  0.4. 
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Fig.  21.  Wind  speed  dependence  of  CKICD  from  this  study  (green  squares)  compared  with 
previous  studies.  Symbols  are  as  in  Fig.  20. 
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The  results  from  this  study  do  provide  some  support  of 
the  similarity  hypothesis  of  Emanuel  (2003)  that  CD  and 
CK  are  constants  at  extreme  wind  speeds.  The  un¬ 
certainty  in  the  coefficients  and  limited  SST  range  pre¬ 
cludes  testing  of  the  hypothesis  that  the  CK  is  SST 
dependent.  These  results  suggest  that  the  lower  bound 
of  the  E-PI  range  for  Hurricane  Isabel  presented  in 
BM08  would  be  the  most  accurate  (cf.  their  Fig.  1). 
These  CKICD  estimates  yield  potential  intensity  esti¬ 
mates  below  50  m  s”1  without  including  dissipative 
heating.  Recent  studies  suggest  that  these  discrepan¬ 
cies  between  observations  and  E-PI  theory  may  be 
largely  due  to  gradient  wind  imbalance  in  the  boundary 
layer  (Smith  et  al.  2008;  Bryan  and  Rotunno  2009).  The 
current  results  are  consistent  also  with  recent  three- 
dimensional  numerical  simulations  that  intensified  to 
major-hurricane  status  with  CKICD  ratios  as  low  as  0.1 
(Montgomery  et  al.  2010). 

These  new  estimates  of  CD  and  CK  may  help  to  pro¬ 
vide  a  basis  to  improve  surface  flux  parameterizations, 
potential  intensity  theory,  and  our  understanding  of 
tropical  cyclone  intensity  change.  Although  there  is  still 
considerable  uncertainty  in  the  magnitude  of  the  air-sea 
fluxes  and  the  corresponding  bulk  exchange  coefficients 
in  tropical  cyclones,  these  estimates  provide  some  new 
evidence  that  CD  and  CK  do  not  increase  at  wind  speeds 
above  50  m  s-1.  Furthermore,  the  current  results  sug¬ 
gest  that  the  CK/CD  ratio  does  not  change  significantly 
from  minor-  to  major-hurricane  wind  speeds.  Additional 
observations  of  air-sea  exchange  at  major-hurricane 
wind  speeds  are  recommended  to  further  reduce  the 
uncertainty  in  these  estimates. 
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APPENDIX  A 

Error  Analysis 

To  test  the  accuracy  of  the  methodology  used  in  this 
study,  momentum  and  energy  budget  analyses  were 
computed  using  idealized  numerical  simulations  from 


the  two-dimensional,  axisymmetric  hurricane  model 
developed  by  Rotunno  and  Emanuel  (1987,  hereafter 
RE87)  and  the  three-dimensional  Advanced  Research 
Weather  and  Forecasting  Model  version  3. 0.1.1  (hereaf¬ 
ter  simply  WRF).  Details  on  the  numerical  configura¬ 
tions  of  the  models  are  given  in  Bell  (2010).  The  initial 
atmospheric  condition  used  for  both  models  was  a 
thermodynamic  background  defined  by  the  Jordan 
(1958)  mean  tropical  Atlantic  sounding  modified  to  be 
in  thermal  wind  balance  (Smith  2006)  with  an  analytic 
mesoscale  vortex  [RE87,  Eq.  (37)].  All  simulations 
were  approximately  steady  state  after  144  h  of  in¬ 
tegration  time. 

The  accuracy  of  the  methodology  was  evaluated  by 
comparing  the  retrieved  exchange  coefficients  from  sim¬ 
ulated  observations  with  known  exchange  coefficients 
derived  from  the  numerical  simulations.  One  WRF 
simulation  and  five  RE87  simulations  were  performed 
with  varying  exchange  coefficients.  A  capped  roughness 
length  of  2.85  X  10  3  m  and  fixed  moisture  roughness 
length  of  1  X  10  6  m  were  used  in  the  WRF  simulation  to 
approximately  match  the  derived  magnitudes  of  CD  and 
CK  from  this  study.  A  variety  of  exchange  coefficient 
values  were  prescribed  explicitly  through  bulk  aero¬ 
dynamic  formulas  in  the  RE87  simulations.  Initial  tests 
indicated  that  the  fluxes  and  exchange  coefficients  could 
be  retrieved  exactly  from  the  simulations  with  all  budget 
terms  included  and  no  observational  errors  (not  shown). 
Estimates  of  errors  in  the  derived  exchange  coefficients 
using  real  observations  were  calculated  by  subsequently 
neglecting  residual  budget  terms  [Eqs.  (4)  and  (8)]  and 
also  by  adding  errors  to  the  simulated  observations.  It 
must  be  noted  that  some  error  was  introduced  by  tem¬ 
poral  averaging  and  interpolating  the  data  from  the 
staggered  to  unstaggered  model  grids,  and  in  the  WRF 
case  by  azimuthally  averaging.  The  interpolation  led  to 
nonnegligible  errors  in  axisymmetric  mass  continuity, 
which  required  correction  by  adjusting  w  through  a  ver¬ 
tical  integration  of  the  mass  continuity  equation  on  the 
temporally  averaged,  unstaggered  grids. 

An  example  of  the  retrieved  CD  and  CK  from  the 
simulated  observations  with  the  residual  budget  terms 
neglected  and  no  observational  errors  is  shown  in 
Fig.  Al.  The  results  are  from  model  output  averaged 
over  eight  “missions”  of  6-h  duration  from  144  to  192  h, 
using  control  volumes  similar  to  those  used  with  the  real 
observations.  One  standard  deviation  from  the  mean  CD 
and  CK  from  the  eight  synthetic  missions  is  indicated  by 
the  vertical  bars.  The  retrieval  methodology  is  able  to 
reasonably  distinguish  between  different  values  of  the 
exchange  coefficients  from  the  six  different  simulations, 
with  differences  from  the  actual  magnitudes  on  the  or¬ 
der  of  50%  or  lower. 
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Fig.  Al.  Retrieved  exchange  coefficients  vs  known  coefficients 
from  synthetic  observations  after  neglecting  budget  residual  terms. 
Gray  circles  indicate  WRF  retrievals,  and  black  squares  indicate 
RE87  retrievals.  Dashed  gray  line  indicates  linear  fit  of  the  data,  with 
the  linear  coefficients  and  correlation  at  the  bottom  of  each  panel. 
Error  bars  indicate  one  standard  deviation  in  the  retrieved  quantity. 


SST  errors  were  assessed  in  a  straightforward  manner 
by  adding  error  to  the  SSTs  from  the  model  simulations 
during  the  CK  retrievals.  Errors  in  the  SST  were  ap¬ 
proximately  linearly  correlated  with  errors  in  the  retrieved 
CK  of  about  20%  °C_1.  We  estimate  the  uncertainty  in  the 
SST  estimates  derived  from  the  satellite  imagery  and 


AXBTs  to  be  on  the  order  of  1°C  or  less,  yielding  a  CK 
uncertainty  of  ±20%  from  this  error  source. 

Center  errors  were  assessed  by  repositioning  the  ob¬ 
servations  and  reprojecting  the  wind  components  in  a 
cylindrical  coordinate  system  displaced  from  the  known 
center.  The  largest  error  was  found  to  result  from  the 
reprojection  of  tangential  wind  to  apparent  radial  wind. 
For  example,  a  5-km  displacement  of  a  60  m  s”1  tan¬ 
gential  wind  at  20-km  radius  can  yield  up  to  a  15  m  s”1 
error  in  the  radial  wind  velocity.  Integrating  a  radial 
velocity  error  weighted  by  mass,  radius,  and  total  energy 
or  absolute  angular  momentum  can  result  in  a  large 
error  in  the  resolved  flux  across  the  boundary  of  the 
control  volume.  Radial  velocity  errors  arising  from  center 
displacement  change  with  radius,  which  also  introduces 
divergence  errors  that  affect  vertical  motion  and  the  re¬ 
solved  flux  through  the  top  of  the  control  volume.  When 
using  only  a  single  radial  pass  for  the  retrieval,  errors  re¬ 
sulting  from  center  displacements  can  yield  up  to  a  ±25% 
error  in  the  retrieved  exchange  coefficients  per  kilometer 
of  displacement.  If  the  center  errors  are  random,  how¬ 
ever,  then  temporal  averaging  will  mitigate  errors  in  the 
analyzed  wind  components.  The  CBLAST  analyses  were 
therefore  averaged  over  multiple  radial  passes  over  ap¬ 
proximately  6-h  periods  to  minimize  this  error.  The 
temporal  averaging  reduces  the  estimated  uncertainty  in 
the  retrieved  CD  and  CK  to  about  20%  for  up  to  3-km 
average  center  errors. 

Analysis  errors  were  assessed  by  comparing  retrievals 
using  SAMURAI  on  the  simulated  observations  with 
small  random  observational  errors.  This  yielded  less  than 
20%  error  in  the  retrieved  fluxes  compared  to  using  the 
interpolated  model  data  directly.  The  retrievals  were  not 
overly  sensitive  to  the  exact  magnitude  of  the  surface 
winds  used  to  convert  the  retrieved  surface  layer  fluxes  to 
bulk  exchange  coefficients  compared  to  other  analysis 
errors.  The  estimated  bulk  Richardson  number  was  near 
zero  because  of  the  strong  vertical  wind  shear  in  the  ana¬ 
lyzed  cases;  therefore,  no  stability  corrections  were  made 
to  explicitly  adjust  the  exchange  coefficients  to  neutral 
stability.  Uncertainties  in  the  surface  wind  derived  from 
the  dropsonde  and  SFMR  data  and  departures  from  neu¬ 
tral  stability  are  therefore  included  in  the  20%  estimate. 

Errors  in  the  resolved  fluxes  are  compounded  by  in¬ 
tegration  over  the  depth  and  width  of  the  budget  control 
volume.  This  suggests  that  thinner,  shallower  control 
volumes  would  be  prone  to  less  numerical  error,  but 
physical  constraints  and  the  data  distribution  must  be 
considered  also.  Optimal  control  volumes  were  some¬ 
what  different  using  the  RE87  and  WRF  simulated  ob¬ 
servations.  The  WRF  simulations  suggested  an  optimal 
depth  of  400  m,  but  the  RE87  simulations  suggested 
a  deeper  volume.  Sensitivity  tests  with  the  simulated 
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observations  indicated  that  aspect  ratios  of  20-30  (width 
divided  by  height)  yielded  the  least  error  in  the  retrieved 
CD  and  CK  from  both  the  RE87  and  WRF  simulations. 
Because  of  uncertainties  in  the  simulations,  the  mini¬ 
mum  depth  of  the  control  volume  was  largely  determined 
by  our  expectation  that  the  unresolved  vertical  eddy 
fluxes  are  small  near  the  top  of  the  boundary  layer.  Wider 
control  volumes  yield  a  larger  integrated  surface  flux 
that  is  more  tolerant  to  errors  in  the  neglected  residual 
terms,  but  are  limited  in  maximum  width  by  the  data 
distribution.  Therefore,  a  compromise  was  struck  be¬ 
tween  physical  and  numerical  constraints  on  the  size 
of  the  control  volume.  A  set  of  72  control  volumes  was 
constructed  for  each  mission  to  maintain  an  optimal 
aspect  ratio  while  also  preserving  the  physical  and  data 
distribution  constraints. 


APPENDIX  B 

SAMURAI  Technique 

The  analysis  uses  an  incremental  form  of  the  variational 
cost  function  that  avoids  the  inversion  of  the  background 
error  covariance  matrix  by  using  a  control  variable  x, 
similar  to  the  forms  in  Barker  et  al.  [2004,  Eq.  (2)]  and  Gao 
et  al.  [2004,  Eq.  (7)]: 

/(x)  —  ixTx  +  ^(HCx  —  d)TR-1(HCx  —  d),  (Bl) 

where  H  denotes  the  linearized  observation  operator,  C 
denotes  the  square  root  of  the  background  error  co- 
variance  matrix,  R  denotes  the  observation  error  co- 
variance  matrix,  and  d  =  y  -  h(xb)  denotes  the  difference 
between  the  observations  y  and  the  nonlinear  observa¬ 
tion  operator  applied  to  the  background  state  estimate 
h(xb).  In  the  current  study,  h  and  H  are  equivalent.  The 
cost  function  is  minimized  using  a  conjugate  gradient 
algorithm  (Polak  1971;  Press  et  al.  2002)  to  find  the  at¬ 
mospheric  state  where  the  gradient  with  respect  to  x  is  0: 

W(x)  =  (I  +  CTHTR_1HC)x-  CTHTR_1d.  (B2) 

We  can  express  the  transform  from  the  control  vari¬ 
able  to  an  analysis  increment  as  an  operator  sequence 
8x  =  Cx  =  SDFx.  The  SDF  matrix  transforms  represent 
the  cubic  B -spline  transform,  standard  deviation  of  the 
background  error,  and  recursive  filter  operators,  respec¬ 
tively.  For  brevity,  only  the  one-dimensional  transform 
for  the  radial  direction  is  illustrated  here.  An  additional 
transformation  in  the  z  direction  follows.  The  spline 
transform  S  is  given  by 


S  =  (P  +  Q)-1 , 

P  =  bmm'E  Pmm'  =  «Mr> A'W  dr, 

3D 

Q  =  [WE  W  =  |fl£,(W»'W  dr ’ 

4>m(r)  =  ^  (  a  ^  ’  ^°r  m  ^  ^  and  r^D, 

(B3) 

where  O  is  the  cubic  B-spline  given  by 

|(2  -  |^  |  )3  -  |(1  -  |^  |  )3  if  1>  |^|  >0, 

if  2>|||>1, 

0  if  HI  >2.  (B4) 

Note  that  P  becomes  the  identity  matrix  for  orthog¬ 
onal  basis  functions,  but  is  required  for  computing  the 
cubic  B-spline  coefficients.  Here  Q  is  a  third  derivative 
constraint  to  reduce  Gibb’s  oscillations  in  the  spline 
transform  (Ooyama  2002),  with  eq  being  a  specified  ef¬ 
fective  spatial  filter  cutoff  length.  The  resulting  in¬ 
crement  or  analysis  can  then  be  evaluated  at  any  point  in 
the  physical  domain  through  the  inner  product  of  the 
basis  functions  and  the  spline  coefficients.  For  example, 
for  the  radial  wind 

u{r )  =  (f)Tx  +  4>Tx,  (B5) 

where  the  second  term  on  the  right-hand  side  represents 
imposed  boundary  conditions  (Ooyama  2002). 

For  the  current  study,  the  background  error  correla¬ 
tions  were  assumed  to  be  Gaussian  and  isotropic,  and 
were  calculated  using  an  efficient  recursive  filter  oper¬ 
ator  that  replicates  the  effects  of  this  correlation  (Purser 
et  al.  2003).  The  operator  combination  DF  is  the  appli¬ 
cation  of  the  background  error  covariance  matrix,  where 
D  is  the  standard  deviation  of  the  background  errors  and 
F  is  the  recursive  filter.  The  standard  deviation  of  the 
background  errors  was  purposefully  set  very  high  given 
no  prior  knowledge  of  background  state  other  than  the 
objective  analysis  from  BM08.  A  large  background  er¬ 
ror  standard  deviation  has  the  detrimental  side  effect  of 
making  the  spline  analysis  unconstrained  in  data-poor 
regions.  The  recursive  filter  length  scale  acts  both  as  an 
effective  distance  for  the  influence  of  the  observations 
and  as  a  spatial  filter.  A  large  length  scale  helps  spread 
the  information  provided  by  an  observation  across  data 
gaps  but  also  removes  finescale  detail  from  the  analysis. 
Sensitivity  tests  indicated  that  a  6A  background  error 
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length  scale  (where  A  is  the  horizontal  or  vertical  grid 
scale)  was  a  good  compromise  between  smoothing  and 
data  density  constraints,  and  was  adequate  for  the  cur¬ 
rent  study. 

The  control  variable  state  vector  is  constructed  in 
axisymmetric  cylindrical  space  and  is  given  by  x  = 
{ prv ,  i//,  S' ,  q Pa}T,  where  p  denotes  the  moist  density 
including  water  vapor,  i/j  denotes  the  transverse  stream- 
function,  S  denotes  the  moist  static  energy  (cpT  +  Lq  + 
gz ),  pa  denotes  the  dry  air  density,  primes  represent  de¬ 
partures  from  a  static  background  reference  state,  and 
the  remaining  symbols  are  as  defined  in  section  2a.  The 
mean  tropical  sounding  from  Jordan  (1958)  was  used  as 
the  reference  state  in  this  study.  Radial  and  vertical  mo¬ 
mentum  were  recovered  from  the  streamf unction  using 
pu  =  —  difr/dz  and  pw  =  di/j/dr.  The  vertical  and  radial 
resolutions  of  the  analyses  used  in  the  study  were  100  m 
and  1  km,  respectively. 
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